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Exciton condensation driving the periodic lattice distortion of 1 T-TiSe 2 
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We address the lattice instability of 1 T-TiSe2 in the framework of the exciton condensate phase. 
We show that, at low temperature, condensed excitons influence the lattice through electron-phonon 
interaction. It is found that at zero temperature, in the exciton condensate phase of lT-TiSe2, 
this exciton condensate exerts a force on the lattice generating ionic displacements comparable in 
amplitude to what is measured in experiment. This is thus the first quantitative estimation of the 
amplitude of the periodic lattice distortion observed in lT-TiSe2 as a consequence of the exciton 
condensate phase. 

PACS numbers: 71.35.Lk,71.45.Lr,71.38.-k 
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In a semimetallic or semiconducting system exhibit- 
ing a small electronic band overlap or gap, the Coulomb 
interaction, when poorly screened, leads to the forma- 
tion of bound states of holes and electrons, called exci- 
tons. If their binding energy Eb is larger than the gap, 
they may spontaneously (without any optical pumping) 
condense at low temperature and drive the system into 
a new ground state with exotic properties. This new 
ground state, called the excitonic insulator phase, has 
been theoretically predicted in the 1960s Since that 
time, among the numerous systems proposed for the re- 
alization of this peculiar phase, only a few turn out to 
be serious candidates [2-6] . Recently we investigated the 
charge density wave (CDW) system 1 T-TiSe2 with angle- 
resolved photoemission spectroscopy (ARPES), favour- 
ing the excitonic insulator phase scenario as the origin of 
the CDW phase @, Q , as suggested earlier 0, [To| . Fur- 
thermore, a superconducting phase has been discovered 
in this material upon copper intercalation [11[ and pres- 
sure 12]. This produces a very interesting phase diagram 
reminiscent of the one of the iron pnictides, in the sense 
that a density wave phase gives way to a superconduct- 
ing dome upon chemical intercalation. The nature of 
the competition between the ordered (CDW) phase and 
the superconducting dome is of central interest. Indeed 
Sawatzky et ai argued that the large electronic polar- 
izabilites of pnictides and chalcogenides, which can be 
interpreted as the consequence of virtual excitonic ex- 
citations, may be at the origin of the superconducting 
phase of this new class of materials, reminding the model 
of Little developped for organics 13l \l4\. 

The quasi two-dimensional material 1 T-TiSe2 under- 
goes a phase transition towards a commensurate 2x2x2 
CDW phase below the critical temperature T c c± 200K. A 
weak periodic lattice distortion (PLD) accompanying the 
CDW (which requires only electronic degrees of freedom) 




FIG. 1: Atomic structure of lT-TiSe2, where the Se (Sel 
and Se2) atoms are in (not exactly regular) octahedral coor- 
dination around the Ti atoms. The orthonormal axis system 
points from the center Ti atom towards the neighbouring Se 
atoms (with a slight deviation of 1.6°). The crystallographic 
axes a.b.c are also shown. 



has been measured, involving small ionic displacements 
< 0.1 A [15]. The occurence of this PLD lead Hughes 
to suggest a band Jahn- Teller effect as the driving force 
of the CDW in lT-TiSe 2 [16]. In this respect, Motizuki 
and coworkers, based on a tight-binding (TB) fit to the 
band structure calculated by Zunger and Freeman [l7| 
found that, by optimizing electronic vs elastic energy, 
the observed CDW is realized for an ionic displacement 
very close to the measured one [18[- However, the small 
ionic displacements in comparison with the high spectral 
weight carried by the backfolded bands as observed by 
ARPES supports rather an electronic origin of the CDW 

0. 

In this perspective, it is crucial to know whether ionic 
displacements of a reasonable amplitude may appear at 
all as a consequence of exciton condensation in the low 
temperature phase. Here, we address this question and 
study the influence of an exciton condensate on the lat- 
tice. 

First, we derive the electron-phonon coupling in the 
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framework of the TB formalism. Focusing particularly 
on the valence and conduction electrons, we derive a for- 
mula relating the ionic displacements to the presence of 
an exciton condensate, the amplitude of which is directly 
related to the order parameter characterizing the low 
temperature phase. Applying this formula to the case of 
1 T-TiSe2 , we calculate the amplitude of ionic displace- 
ments. We find values similar to those obtained from 
experiment. This demonstrates that the exciton conden- 
sate phase, as a possible origin of the CDW phase of 
lT-TiSe2, can also account for the PLD. 

The TB formalism for the electronic band structure 
and the coupling to the lattice described below is similar 
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to that developed by Yoshida and Motizuki 
is applied to the structure of lT-TiSe2, which consists 
of planes of Ti atoms forming a triangular lattice. Each 
of these Ti atoms is in octahedral coordination with its 
six neighbouring Se atoms. Then the crystal consists of 
a regular stacking of such Se-Ti-Se layers along the c 
direction. For the TB calculations of the present study, 
a cluster of atoms centered around one Ti atom will be 
considered. The orthonormal axis system is as shown in 
Fig. [U In our calculations, for the Ti atoms, we include 
the five 3d-orbitals of xy, yz, xz, x 2 — y 2 and 3z 2 — r 2 
symmetry and the three 4p-orbitals for each of the two 
Se atoms (Sel and Se2), of x, y and z symmetry. In total 
we have 11 orbitals. The TB electronic Hamiltonian then 
reads 

Hei = J2 J2 J ^ ~ ^ +^-t v ) 

IV pv a(3 
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Here Ri,Ri> are vectors of the Bravais lattice and r a ^r v 
are the positions of the ions (Ti, Sel or Se2) labelled 
/i, v inside the unit cell. The indices a, j3 label the 11 
orbitals and n is the index of the bands in which the 
operators create electrons. The transfer matrix J, 
with the eigenvectors T, consists here only of two-center 
integrals for simplicity. 

We now introduce ionic displacements of the form 



the electron-phonon interaction 
Hei-ph = ^2 ^29nn'(k,q, \)c ] n {k)c n > {k - q)Q(q, A) (3) 

nn ' kq,X 

where the electron-phonon coupling constant 

p pv a(3 

xW(l-g>-^. [e(/i,£A)e 



dx 



X=p+T^-T u 

e>,<?,A)], (4) 



with p = Ri—Ri>, involves the derivatives of the transfer 
matrix dJ a p / dx. 

In what follows, we focus on the influence of excitons. 
In our simplified model of the band structure of 1 T- 
TiSe2, we essentially consider the topmost valence band 
having its maximum at the center of the Brillouin zone 
(r point) and the three symmetry equivalent conduction 
bands at the border of the Brillouin zone (L points), 
whose extrema are separated from T by the wave vec- 
tors Wi (i = 1,2,3) j8|. These excitons are composed of 
holes created by a(k) in the valence band (near its max- 
imum) with a wave vector k and electrons created by 
b\(k) in the conduction band i (near its minimum) with 
wave vector k + Wi. Thus the sum over the band indices 
n, n' in equation [3] is restricted to terms mixing a and b 
operators only, so that 

Hei-ph = ^ ^ Q(q,\)g abi (k,q,\)a\k)bi(k-Wi - q) 

1 kq,X 

+Q{& X)g hia (k, q, X)b\(k - Wi)a(k - q) 
= Yl X )dabi (p, ~Wi, \)a ] (p)bi(p) 

i p,X 

+Q(wi, X)g hi0b (p + Wi,Wi, X)b\(p)a(p). 

We considered only q = — Wi in the term involving g a \ )i 
and q = Wi in that involving gb ia (together with the sub- 
stitution p = k — Wi), restricting ourselves to the scatter- 
ing between the extrema of the bands. Then, averaging 
(H e i- P h)ei to lowest order over the electronic degrees of 
freedom yields the contribution of condensed excitons to 
the phonon Hamilonian 
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Q(-Wi, X)g a bi (p, -Wi,X) (a f (p) bi (p) ) 



where M M is the mass of the ion labelled /i, e a polariza- 
tion vector and Q the normal coordinate of the phonons. 
Here (<f, A) is the ionic displacement for the atom la- 
belled ji (Ti, Sel or Se2) and associated to a particular 
mode <f, A. Equation [2] provides us with a direct way to 
compute the amplitude of the displacement of each ion, 
once we get a value for Q. This is our goal in the next 
paragraphs. After introducing the small ionic displace- 
ments in the argument of the transfer matrix J in 
equation [TJ we expand J to first order in ui^ to deduce 



+Q(wi, X)g bia (p + Wi,Wi, \){b\(p)a(p)). 



In analogy to the BCS-theory, the averages (b\a) are 
related to anomalous Green's functions Fi{p,r) = 
(—i)(Tbl(p,r)a(p)) (introduced in reference [8]), so that 
the previous equation becomes 

(H e i_ ph ) e i = «5^^Q(-^i, X)g abi (p, -Wi, \)f}(p,0)} 

i p,X 

+Q(wu X)gb z a(p + m,Wi, X)Fi(p, 0) 

= Hph — x- 
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This exciton-phonon Hamiltonian H p h- X can be further 
simplified using the inversion symmetry of the system to 
replace —wi by Wi and using also the property Fi(p, 0) = 
F-(p,0), giving 

Hph-x = i ^ Q(w, A) Fj(p, 0) 

x [^(p,^, A) +^a(p + ^,^i,A)] . 

From this last equation, the equilibrium condition for 
the lattice in the presence of a condensate of excitons, 
d(H p h,o + Hph-x) /dQ(wi, A) = 0, leads to an expression 
for the normal coordinate of the phonons Q caused by 
the exciton condensate 



Q{wi,\) 



(5) 



where H ph ,o = (1/2) £\ A u 2 (wi, X)Q*(wi, \)Q(wi, A) is 
the bare Hamiltonian of the lattice (in the absence of 
the exciton condensate). By using equation [2] we can 
relate Q to the resulting ionic displacements. 

Having now an analytical formula for the ionic dis- 
placements thanks to equations [2j 2J [5j we can look for 
the necessary numerical parameters for the final compu- 
tation. We start with the transfer matrix J a p. According 
to Slater and Koster [20[ , its elements are computed as a 
combination of direction cosines (defining the direction of 
the bonds joining two atoms) and transfer integrals (spe- 
cific to the orbit als involved). In our case, these trans- 
fer integrals are determined by fitting a band structure 
computed with density functional theory (DFT). This 
first-principles band structure has been calculated using 
the full potential augmented plane wave plus local or- 
bitals (APW+lo) method with the generalized gradient 
approximation in the parametrization of Perdew, Burke 
and Ernzerhof [21], in the local density approximation, as 
implemented in the WIEN2K software package [22[ . The 
numerical and unit cell parameters [23L l24j are similar to 
those used in reference [25j. 

The resulting DFT band structure is plotted in Fig. 
El The goal of this DFT calculation (such a DFT band 
structure is discussed in detail in references [25|, |26|) is 
to offer a basis for the TB fit we perform for determining 
the transfer integrals. A least-square fit method was used 
on our DFT band structure along the high symmetry 
directions YMKTA and the obtained transfer integrals 
are summarized in Tab. U j^- The TB fit of the DFT 
band structure is good for states below the Fermi energy 
Ep but becomes poor for states above Ep. This TB 
parametrization allows to build an effective Hamiltonian 
which can be diagonalized (at each k point) to provide 
the eigenvectors T a/Xjn (fc). The derivatives of the trans- 
fer matrix dJ a pjdx, appearing in the electron-phonon 
coupling g nn '(k, q, A), lead to derivatives of the direction 
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FIG. 2: Comparison of the band structure of lT-TiSe2 calcu- 
lated with density functional theory and its fit within a TB 
approach. 



cosines and derivatives of the transfer integrals. The lat- 
ter are evaluated as in Ref. [18j . 

According to the experimental result of Di Salvo et ai 
[lH, we fix the phonon polarization vectors e(/i, <f, A) in- 
volved in the CDW to the direction perpendicular to their 
respective q vector, lying in the ab basal plane. This way, 
only the transverse phonon mode X tr will be considered 
in the following calculations. 



TABLE I: Transfer integrals (in eV) for lT-TiSe 2 , obtained 
from TB fits to first-principles band structure £de^Sd 1 are 
the orbital energies) [291 ]. 

Transfer integral Energy Transfer integral Energy 



t(ppa) 


0.77 


t(pdn) 


0.70 


t(ppn) 


-0.054 


t(ppa) 2 


0.63 


t(dda) 


-0.35 


t(pp7r)2 


-0.028 


t(dd7v) 


0.074 


t(ppa) 3 


0.61 


t(ddS) 


-0.049 


t(pp7r) 3 


-0.096 


t(pda) 


1.3 






e p = -2.0 eV 


£d£ = 0.74 eV 


e dl = 1.2 eV 





The anomalous Green's function r = 0) appear- 
ing in equation [5] is calculated as the Fourier transform 
of Fi(jp,z) given by (see reference |8|) 



A(z - 4 +1 (P + w 2 ))(z - 4+ 2 (p + w 3 )) 



V{p,z) 

with the denominator being 

V(p,z) = (z-e v (p))Y[(z-e l c (p-\-Wi)) 

i 

i j^i 

Here, the order parameter A describes the intensity of the 
exciton condensate in the low temperature phase. The 




100 150 200 250 
Temperature (K) 



FIG. 3: . Amplitude of the Ti ionic displacement for a single 
Wi, UTi(wi, Xtr) (the diamonds show the calculated values). 



functions e v and e % c describe the dispersions of the valence 
band and of the three symmetry equivalent conduction 
bands (i = 1,2,3) close to their maximum situated at T 
and at L, respectively. The anomalous Green's function 
Fi is sensitive to the energies appearing in these disper- 
sions near their extrema. We therefore cannot use the 
TB dispersions, which are too rough with this respect 
(however they are essential for the more global treatment 
needed to obtain the transfer integrals), but we need the 
formulas for e v and e % c determined in a previous study 
from fits to ARPES data Q. 

Finally, combining equations [2] and [5j the amplitude of 
the ionic displacement for a single Wi (and for transverse 
phonons) gets the following form 
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EifW>0) [ga bi +g bi a] 



(6) 



We now focus to the particular case of the Ti atoms, 
so that Mti describes the Ti atom mass. For the order 
parameter appearing in the anomalous Green's function, 
we consider a mean-field like temperature dependence 
of the form A(T) = A ^l - (T/T c ) 2 with T c = 200K, 
where the zero value Ao = 115 meV has been determined 
in our recent temperature dependent ARPES study [lo| . 
In formula [6j we use a value of oj{w^ Xtr) =6.3 THz. In 
fact, it is the value estimated by Holt et at. at T ^ 150K 
for the transverse phonon mode, which softens at L at 
the transition [31]. This corresponds to a situation where 
the ionic displacement ui^ is not too large, such that our 
first order development of the electron-phonon coupling 
remains valid. Furthermore, the region close to T c is 
avoided, where anharmonicities cannot be neglected in 
the bare Hamiltonian for the lattice H p h,o- 

Now all the necessary quantities to compute the ampli- 
tude of the ionic displacements in equation [6] are known. 
Fig. [3] summarizes our numerical results. It shows a clear 
temperature dependence, following closely the behaviour 
of the order parameter. Extrapolated to the lowest tem- 
perature, it reaches the value of u^f°(wi, X tr ) = 0.025 
A. Di Salvo et al. inferred from neutron diffraction ex- 
periments a displacement (also for a single-^) of about 



u e *v(wi, X tr ) = 0.04 A at 77K [15]. Therefore our value, 
although being about 60% of the experimental one, re- 
produces the measured ionic displacement for Ti atoms 
within the correct order of magnitude, which is a sub- 
stantial result, considering the approximations made in 
this calculation. The uncertainty on uj(wi,Xt r ) used in 
equation [6] may enhance or reduce this value by a fac- 
tor of 2-3, but it still remains within the correct order of 
magnitude, in agreement with the main message of this 
letter. 

In conclusion, we addressed the question of the appear- 
ance of a periodic lattice distortion in lT-TiSe2- Previ- 
ously we gave strong support for the exciton condensa- 
tion as a purely electronic mechanism responsible for the 
CDW phase in this material [7|, |8[. In this work, we 
elaborate in a tight-binding formalism a formula for esti- 
mating the ionic displacements produced by the presence 
of this exciton condensate through the electron-phonon 
coupling. The calculated amplitude of these ionic dis- 
placements is, at low temperature, of the same order of 
magnitude as what is experimentally found. This is thus 
the first quantitative estimation of the amplitude of the 
PLD observed in 1 T-TiSe2 as a consequence of the ex- 
citon condensate phase. More generally, this result de- 
scribes quantitatively how an excitonic insulator phase 
can give rise to a PLD through electron-phonon interac- 
tion. 
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